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Abstract. The diquark and nucleon are studied in a non-local NJL model. We solve the relativistic 
Faddeev equation and compare the results with the ordinary NJL model. Although the model is 
quark confining, it is not diquark confining in the rainbow-ladder approximation. We show that the 
off-shell contribution to the diquark T matrix is crucial for the structure of the nucleon: without its 
inclusion the attraction in the scalar channel is too weak to form a three-body bound state. 



INTRODUCTION 

The NJL model is a successful low-energy phenomenological model inspired by QCD 
[1]. It has also been applied to the nucleon, see Refs. [2]. However, the lack of con- 
finement makes the model questionable in the baryonic sector. It has been shown that a 
non-local covariant extension of this model can lead to quark confinement for accept- 
able values of the parameters [3]. This occurs due to the fact that the quark propagator 
has no real poles, and consequently quarks do not appear as asymptotic states. There 
are several other advantages of the non-local version over the local NJL model: the non- 
locality regularises the model in a manner such that anomalies and gauge invariance are 
preserved and the momentum-dependent regulator makes the theory finite to all orders 
in the l/N c expansion. Finally the dynamical quark mass is momentum dependent in 
contrast to the ordinary NJL model and consistent with lattice simulations of QCD. As 
a result, the non-local version of the NJL model may have more predictive power. 

Here we shall use the covariant quark-diquark formalism in a relativistic Faddeev 
equation for the nucleon, where we include only scalar diquark correlations. Due to 
the separability of the non-local interaction, the Faddeev equations can be reduced to a 
set of effective Bethe-Salpeter equations. This makes it possible to adopt the numerical 
method developed for such problems by Oettel et al. [4]. 



THE MODEL 

We consider a non-local NJL model Lagrangian with SU (2) f x SU (3) c symmetry. There 
exists several versions of such non-local NJL models. Regardless of what version is 
chosen, a Fierz transformation allows one to rewrite the interaction in either the qq or 
qq channels. Here we truncate to the scalar (0 + ,T = 0) and pseudoscalar (0~,r = 1) 



mesonic channels. To investigate the nucleon in the quark-diquark picture, one also 
needs to know the qq interaction. We truncate to the scalar (0 + ,T = 0) colour 3 qq 
channel, 



J s (x) = / d 4 X2d 4 X4f(x2-x)f(x-X4)y/ T (x 2 )[C l y 5 T 2 p A ]Y(x4)- (1) 



Here T a = (l,?'y 5 T) and m c is the current quark mass of the u and d quarks. The ma- 
trices /3 A = a/3/2A a (A = 2,5,7) project onto the colour 3 channel with normalisation 
tr(/3 A /3 A ') = 38^' and the t,-'s are flavour SU(2) matrices with tr(r ; T 7 ) = 2% The ob- 
ject C = ifzYs is the charge conjugation matrix. Since we do not restrict ourselves to 
specific choice of interaction, we shall treat the couplings g s and g K as independent. For 
simplicity, we assume the form factor f(x — Xi) to be local in momentum space, since it 
leads to a separable interaction in momentum space. 

The dressed quark propagator S(k) is now constructed by means of a Schwinger- 
Dyson equation (SDE) in rainbow-ladder approximation. Thus the dynamical con- 
stituent quark mass, arising from spontaneously broken chiral symmetry, is obtained 
as [the symbol Tr denotes a trace over flavour, colour and Dirac indices] 



Following Ref. [3], we choose the form factor to be Gaussian in Euclidean space, 
I(pe) = ex P(~ pI/A 2 ), where A is a cutoff of the theory. If one assumes that A is 
related to the average inverse size of instantons 1/p, its choice thus parametrises non- 
perturbative properties of the QCD vacuum. The choice (2) respects Poincare invariance 
and leads to quark, but not colour, confinement, when the dressed quark propagator has 
no poles at real p 2 in Minkowski space. This occurs for 



For large enough values of the dynamical quark mass, the quark propagator has no real 
poles and quarks do not appear as asymptotic states. The propagator still has infinitely 
many pairs of complex poles, both for confining and non-confining parameter sets. This 
is a feature of such models and due care should be taken in handling such poles, which 
can not be associated with asymptotic states if the theory is to satisfy unitary. 

Our model has four free parameters: the current quark mass m c , the cutoff (A), the 
coupling constants g n and g s . We fix the first three to give a pion mass of 136.6 MeV 
with decay constant of 92.4 MeV, as the value of the zero-momentum quark mass in 






(2) 




(3) 



TABLE 1. The model parameters for dif- 
ferent sets, fitted as discussed in text. The 
resulting value of the dynamical quark mass 
m(0) also is shown. 

Parameter set A set B set C 

mo(0)(MeV) 250 350 400 

m(0)(MeV) 297.9 406.2 461.3 

m c (MeV) 7.9 14.4 17.6 

A(MeV) 1046.8 723.4 638.1 

^(GeV- 2 ) 31.6 89.0 132.6 



chiral limit rao(0). We analyse three sets of parameters, as indicated in table 1. Set A is 
non-confining and sets B and C are confining (i.e., they satisfy the condition Eq. (3)). 
The parameter g s is not constrained by these conditions, which allows us to analyse the 
coupling constant dependence through the ratio r s = g s /gn- 



DIQUARK AND NUCLEON SOLUTION 

In the rainbow-ladder approximation the diquark T -matrix can be written as 

T(pi,P2,P3,PA) = f(Pl)f(P3)[r 5 CT2P A ]t(q)[C- l y 5 T 2 PV(P2)f(P4) 

X8(pi+P 2 -P3-P4), (4) 

where q = p\— p-i = Pa~ Pi and 

t{q) = y (5) 

Js{q) = /Tr|^/ 2 (-^)[ r5 CT 2j 8 A ]5(-^) r [C-V5T2j8l% + ^)/ 2 (? + ^).(6) 

In the above equation the quark propagators are the solution of the rainbow SDE Eq. (2). 
The diquark bound state is located at the pole of the T-matrix. The denominator of 
Eq. (5) is the same as that in pionic channel, thus one may conclude that at r s = 1 
the diquark and pion are degenerate. This puts an upper limit to the choice of r s , since 
diquarks should not to condense in vacuum. 

The relativistic Faddeev equation can be written as an effective two-body BSE be- 
tween a quark and a diquark due to the separability of the two-body interaction in 
momentum-space. Using the the on-shell approximation this becomes an eigenvalue 



problem [4], 



/ ^4 K ar(P^P)v r p(KP)=X(P 2 )^ a p(p,P) 1 (7) 
il/(p,P)=S(p q )D(p d )<l>(p,P), (8) 
K(p,k;P) = -^-X{P x )S T {q)X{P 2 ). (9) 

&dqq 

Here ^(p^P) is the nucleon bound-state wave function related to the vertex function 
(j)(p,P) by truncation of the legs. The parameters p and P are the relative and total 
momenta in the quark-diquark system, respectively. We define p q = r\P + p, p d = 
(1 — r\)P — p and q = —p — k + (1 — 2r\)P. The relative momentum of quarks in the 
diquark vertices X and X are defined as p\ = p + k/2 — (1 — 3r\)P/2 and p 2 = —k — 
p/2 + (l — 3 r/) P/2 respectively. (For r] see below.) 

These equations should be solved for the largest eigenvalue of k(P 2 ) with the con- 
straint that A(P 2 ) = l/g 2 dqq at P 2 = Mjy, where Mjv is nucleon mass. Here g dqq denotes 

the diquark-quark interaction coupling, g d 2 q = (d k iJ s ) \ k 2 =M i, where J s is defined in 
Eq. (6). In the usual on-shell approximation the scalar diquark propagator is taken to 
be D(pd)~ l = —(p d + Mj) where M d is the pole of the scalar diquark T-matrix. The 
functions X and its adjoint X = JoX^Jq are the vertex functions of two quarks within the 
diquark and can be obtained from Eq. (4), 

X( Ph k d ) =g dqq Y5CT 2 P A f(pi + (t-o)k d )f(-pi + ok d ). (10) 

The non-locality of our model naturally leads to an extended diquark, and provides a 
sufficient regularisation of the ultraviolet divergences for a diquark-quark loop. Above 
we have introduced the Mandelstam parameters r\ and o, T7, c e [0,1]. In principle, the 
eigenvalues should not depend on these parameters if the formulation is Lorentz covari- 
ant. However, in practise due to numerical errors and singularities in the propagators, 
one usually only finds that there are substantial plateaus where the results are r\- and 
<7-independent. The nucleon vertex function <j>(p,P) should correspond to a state of 
positive energy, positive parity and spin- 1/2. This leads to the general structure 

<j,(p,P) = ( Si + ] ^p'S 2 ^A\ (11) 

where = 5(1 + jj^-) is a positive energy projection for a nucleon with mass M#, 
and S\ and S2 are scalar function. We now solve the equations in the rest frame of the 
nucleon. We expand the scalar functions in terms of Chebyshev polynomials of the first 
kind, which turns out to be very efficient for such problems [4]. 



RESULTS 



The scalar diquark mass Md is obtained as a pole of Eq. (5). We find that for a wide 
range of r s , regardless of the parameter set and confinement, a bound scalar diquark 




FIGURE 1. (a) Shows the diquark masses with respect to r s for three sets of parameters, (b) Shows the 
discrepancy between the on- and off-shell approximation for three diquark masses (parameter set B). In 
(c) the nucleon mass is shown as a function of diquark mass for set B, the values are given in GeV. 



exists. The diquark masses for various value of r s for different parameter sets are plotted 
in Fig. la. For sets A, B, C no bound scalar diquark exists for r s < 0.110,0.112,0.146, 
respectively. Please note that the confinement of the diquark for very small r s is due to the 
screening effect of the ultraviolet cutoff and should not be associated to the underlying 
QCD confinement which stems from infrared divergence of gluon and ghost propagators. 
Having said that it is possible that diquark confinement may arise beyond the rainbow- 
ladder approximation [5]. 

In order to solve the Faddeev equation numerically, we have first used the on-shell 
approximation. We observe that one can not generate a three-body bound state in this 
model, without an artificial enhancement of the quark-diquark coupling of about 3. A 
similar feature has been observed in the ordinary NJL model as well [2], but the effect 
is even more severe here. As can be seen from Fig. 1 a decrease in r s leads to a larger 
diquark mass, and an increase in the off-shell contribution to the qq T-matrix. It is this 
off-shell contribution that leads to a bound nucleon, and they are indeed substantial as 
can be seen in Fig. lb. A preliminary result is shown in Fig. lc. We also show a fictitious 
diquark-quark threshold defined as M qq + m q , where m q is the real part of the first pole 
of quark propagator. Given this definition the diquarks in the nucleon are much more 
loosely bound in the non-local model than in the ordinary NJL model. Nonetheless, we 
obtain a bound nucleon solution near its experimental value. 
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